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Abstract 

We construct an asymptotic series for a general solution of the Einstein equa- 
tions near a sudden singularity. The solution is quasi isotropic and contains nine 
independent arbitrary functions of the space coordinates as required by the struc- 
ture of the initial value problem. 

PACS number: 98.80.-k 



1 Introduction 

The standard homogeneous and isotropic Friedmann-Lemaitre cosmological models of 
general relativity are determined by the solution of two independent Einstein equations 
for three unknown time-dependent functions, the metric scale factor a(t), the fluid den- 
sity, pit), or the isotropic fluid pressure, pit). If an equation of state p = f(p) is chosen, 
the system closes and the two remaining unknown functions are determined. Physically 
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reasonable equations of state, like those for perfect fluids, produce well-behaved expan- 
sion factors with behaviours that offer simple Newtonian interpretations. However, the 
challenge of providing a compelling explanation for the observed acceleration of the uni- 
verse in terms of a 'dark energy' fluid has led to an exploration of other less familiar 
equations of state [JJ, motivated by the form of bulk viscous stresses [2j[3lll], where dp/ dp 
is not everywhere continuous, or scenarios in which there is no equation of state at all. 
These less constrained material stresses introduce new possibilities for the evolution of 
Friedmann universes and can produce unexpected types of finite-time singularity which 
are far weaker than the curvature singularities that often accompany geodesic incom- 
pleteness [5j EJ [7J [8j [9j [10] . We found [TTJ [12j [13] a whole new class of pressure-driven 
singularities that keep the scale factor, a, the cosmological expansion rate, a /a, and the 
density, p, finite while permitting the pressure, p, and the deceleration —d/a, to blow 
up at a finite time despite the energy conditions p > and p + 3p > holding, and 
regardless of the three-curvature. There is no geodesic incompleteness associated with 
these singularities. 

The stability of general relativistic solutions with sudden singularities in the presence 
of small scalar, vector, and tensor perturbations has been demonstrated in a gauge 
covariant formalism [14], and they have also been shown to be stable against quantum 
particle production processes as the singularities is approached |15J. They have been 
investigated in a number of different cosmological setting and their behaviour has been 
classified in the light of other types of finite-time singularity that can arise in isotropic 
and anisotropic cosmologies rjllltl^tlEltlHJEniE^ The formal 

definitions of these singularities have been discussed by Lake |27] in the case where 
infinities occur in second derivatives of a(t), but similar examples exist where the sudden 
singularity occurs in (arbitrarily) higher derivatives of a(t) and so no energy conditions 
are threatened [TBI ITS] . Recently, we have also provided another formal construction 
procedure for these solutions in isotropic universes using fractional power series [28, 29J. 

Sudden singularities have appeared in a wide range of cosmological investigations 
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and our aim in this paper is to try to obtain information about the general cosmological 
behaviour in the neighbourhood of such a singularity by an asymptotic analysis of the 
Einstein equations without using simplifying assumptions of homogeneity or isotropy. We 
shall use a function- counting method, familiar from ref [30J. Past uses of this method 
have failed to identify general behaviour near isotropic singularities [301 Ell E2] because of 
the ubiquity of chaotic behaviour near curvature singularities. However, although the full 
complement of independently arbitrary functions did not appear, the series expansions 
obtained provided useful asymptotic forms for particular cosmological models containing 
fluids and combinations of fluids [33J EI]- The only use of this technique which identified 
a general behaviour in the Einstein cosmological equations is the study of the late-time 
behaviour of expanding universes with p > and p + 3p > in the presence of a positive 
cosmological constant. A series expansion around the de Sitter solution [35] produces a 
general complement of arbitrary functions in accord with the expectations of the cosmic 
no hair theorem ES E3 ESI E2] • 

In this paper we are going to seek a series expansion around the particular FRW 
solution with a future sudden singularity introduced in ref [HI [12]. In the limit that the 
sudden singularity is approached at finite future time we will determine the number of 
independently arbitrary functions of the three space variables that are necessary and suf- 
ficient to prescribe the solution locally on a spacelike hypersurface of constant comoving 
proper time. 

2 Friedmann-Lemaitre sudden singularities 

Consider the most general form of sudden finite-time singularity of the sort discussed 
in [TT| [T2] that can arise in an Friedmann expanding universe. In [TT] we showed that 
singularities of this sort can arise from a divergence of the pressure, p, at finite time 
despite the scale factor, a(t), the material density, p, and the Hubble expansion rate, 
H = a/a remaining finite. The pressure singularity is accompanied by a divergence in 
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the acceleration of the universe, a, at finite time. Remarkably, these singularities occur 
without violating the strong energy conditions p > and p + 3p > 0. They can even 
prevent a closed Friedmann universe that obeys these energy conditions from attaining an 
expansion maximum [12]. In order to prevent the occurrence of these sudden singularities 
it is necessary to bound the pressure. A sufficient condition is to require dp/ dp to be 
continuous or p/ p to be finite. 

We will construct a quasi-isotropic ansatz for the general behaviour by first consid- 
ering the Friedmann universe with curvature parameter k, zero cosmological constant, 
and units 8nG = c = 1. The Einstein equations reduce to 



3# 2 = p-% (1) 
cr 

p + 3H(p + p) = 0, (2) 



In [TT] we constructed an explicit example by seeking, over the time interval < t < t s 
a solution for the scale factor a(t) of the form 



a{t) = l + Bt q + C{t s - t) n , (4) 

where B > 0, q > 0, C and n > are free constants to be determined. If we fix the zero 
of time by requiring a(0) = 0, so Ct n s = — 1, we have 

a(t) = (1 )*(a. - 1) + 1 - (1 - i (5) 

"8 "8 

where a s = a(t s ). Hence, as t — > t s from below, we have 

a q(q - l)Bt^ - t2 ^l'L ~+ ( 6 ) 

whenever 1 < n < 2 and < q < 1; the solution exists on the interval < t < t s . Hence, 
as t — > t s we have a — >■ a s , H — )■ if s and p — >■ p s > where a s ,H s , and p s are all finite but 
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p s — > oo. The infinite deceleration in fl6]) is why this sometimes also termed a 'big brake' 
singularity. By contrast, as t — > there is an initial all-encompassing strong-curvature 
singularity, with H — > oo, p — > oo and p — > oo. From ([6]) and ([3]), we see that p and 
p + 3p remain positive throughout the evolution but, because p is finite asymptotically, 
the dominant-energy condition, \p\ < p, must always be violated as t — )■ t s , [TTJ [27J 116] - 
Using the isotropic and homogeneous solution Q as a guide we consider a series 
expansion of the metric in the neighbourhood of the sudden singularity in a general 
cosmological solution to the Einstein equations. 

3 General series expansions near the singularity 
3.1 Fractional power series 

We use the Landau-Lifshitz [30] notation: Latin indices for spacetime components, Greek 
for spatial components. The general form of the metric in synchronous coordinates is 

ds 2 = dt 2 — r ) a pdx OL dx 13 . (7) 

On approach to the sudden singularity as t — > t s , the solution ()5]) has the linear asymp- 
totic form 

a s + q(l-a s )(l- — ). (8) 

Guided by this FRW example, (EI), and its asymptotic form, (JED, we carry out a 
change of time coordinate to place the future sudden singularity at t = and then 
consider a general series expansion about it, with 

7a/3 = a a p + b afj t + c afj t n + 

where 1 < n < 2, and a a/ 3, b a/ 3, c a/ 3 are functions of the space coordinates only; the inverse 
metric tensor is 
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so that •ja/31 13 ' 1 = ^ 2t- We note that a^aP 1 = 5% and the indices of b a p, c aj 3 are raised by 
a Q/3 . The extrinsic curvature, its derivatives, and its contractions are, to leading order, 

K% = b% + nc% t 11 - 1 + ■■■ , 

K = b + nc t n_1 + • • • , 

= n(n - l)c^- 2 + • • • , 

8K 

— = n(n-l)cf- 2 + -, 

K^K% = b^ + 2nb^t n - 1 + --- , 

if-Kjg = 66^ + n{cb% + 6c^) t"" 1 + • • • 

The components of the Ricci tensor are 

_ IdK l K Ka _ n(n-l)c 2 1 n n _! 

~~ 2 <% 4 a p ~~ 2 4 ^ 2 M ' 

where Q a/ 3 is the Ricci tensor associated with 7 Q/ g, and P a/ 3 is the Ricci tensor associated 
with a a p. In the synchronous reference frame 

9oo = 1, 9o a = 0. 
The 4-velocity u l = (u°,u a ) has u° = uq and 

1 = UiU 1 = u 2 Q - (a a/3 + b a/3 t + c a/3 t n H )u a u p . 

The stress energy tensor for a perfect fluid 

T* = 8n[(p + p^Uj - p5*\ 
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gives (to lowest order) 



T ° = ( p + p )ul-p^p, 

T a = (P + P) U U a ~ (P + P)u a , 



From the Einstein equations 



we use the (°) component to calculate the energy density and the trace of the (^) 
equations to calculate the pressure. The 4-velocity can be calculated then from the (°) 
component. Also we use the (^) components to get restrictions on the arbitrary functions 
of a at 3, b a /3, c a p. We have 

i6n P = I(i2g - K) = (p + b2 - b ^ b ^ _ VLqrcr - be) r- 1 + • • • , 

while from the trace of the component we get 
The Ricci scalar is 

i? = ijj = _„(„ - l)ct n - 2 - ^ + / + 4P - + 6c) t n ~ l + ■■■, 

the (^) components of the Einstein equations give 

R1= l -{R-l^p)5l 



and therefore 



n(n-l)c xa 2 / 6^f4P\ nbc x 
6 6 P l J ^ 6 + 



From the 0(t n 2 ) terms, we get the constraints 



and so the components c a p are all but one determined from a Q/ g. 
Since 7^ is to order t n , the extrinsic curvature will be 0{t n ~ l ) 



(9) 



and the Ricci com- 



ponents .Rq, Rp will be 0(t n ~ 2 ). Hence this will be the only restriction on a a p, b a p, c a p. 
Therefore, we have 6 + 6 + 1 = 13 independent functions. Subtracting the 4 coordi- 
nate covariances which may still be used to remove four functions, leaves 9 independent 
arbitrary functions of the three space coordinates on a surface of constant t time. 

This is the maximal number of independent arbitrary spatial functions expected in a 
local representation of part of the general solution of Einstein's equations near a sudden 
singularity. In general, we expect there will be 6 x g aP and 6 x g aP , plus 3 free velocity 
components, w a ,and p and p, giving a total of 17 independent functions. We can remove 
four of these by using the R° a constraints and four more by using the general coordinate 
covariances. This leaves a total of 9 free functions expected in the general solution for 
the metric locally. If an equation of state had been assumed to relate the pressure to the 
density this number would have been reduced by 1 to 8. 
We have found that the pressure is 0(t n ~ 2 ), so p + p = 0(t n ~ 2 ), and from 



R a = 8vr(p + p)u u a pa 8n(p + p)u, 



since R° a 



O(t ), we get that u, 



a 



0(t ( n 2 )), and the full expansion is 





The equation for the stress energy conservation T 13 . { = is split into time and spatial 
components as follows: 



T i0 .. = T a0 a + d t T 00 + \{KT m + K aP T aP ) = 0, 
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where 

The time component gives 

[(p + p)u% tt + d tP + l -K{p + p) = 0. (10) 

To lowest order, the first term of eq. f llOp is constant (0(t°)). The other two terms to 
lowest order are 0(t n ~ 2 ); therefore, to lowest order the equation is 0(t n ~ 2 ). Hence the 
0(t n ~ 2 ) term of d t p must cancel the 0(t n ~ 2 ) term of \K(p + p). The first term is 

16tt( P + p)u a = -a^(6j iAt - 6 ;/9 ) - na a ^ ;fl - c^t™" 1 + • • ■ 

The second is 

Ifordtp = - n{n ~ l \ b^ - bc)t n - 2 + ■■■, 

and the third is 

1 -K(16np+16np) = - n{n - 1)bC r~ 2 + .... 
2 3 

Therefore this equation will be satisfied to order 0(t n ~ 2 ) if 

n ±^l {K ^ - be) + = 0, 

that is, by 

b^c?" = -be, 

which is indeed satisfied because of the constraint eq. on c a p found earlier. Therefore, 
the time component of T %3 ^ = does not give any new constraints on the terms of the 
series. 

Similarly, the spatial part of the stress-energy conservation is expressed as 

7 Q/3 P;« + d t [(p + p)u p ] + \k{ 9 + p)u fi + K?(p + p)u a = 0. 

The lowest order of the first two terms is 0(t n ~ 2 ) while the lowest order of the last two 
terms is 0(t°) and so they are neglected. For the first term, 

ltorf*p ia = - 2<n ~ l) a^c, a r- 2 + • • • , 
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while for the second, 

16nd t [(p + p)uP] = -n(n - l)a^(c^ - c, u )t n ~ 2 + ■■■ . 
Hence, we must have 

l) a ap c, a + n(n - l)a^(c^ - c. v ) = 

or 

which again is satisfied by the restriction of eq. fl9]) on c a p. Hence, both the temporal 
and the spatial components of the conservation of the stress energy tensor are satisfied 
to lowest order by the constraint on that we have found already and there remain 9 
independent functions in our series solution. 
As an aside we note that he speed of sound is 

v 2 = ^P = P = 4(2 - n)c . 
Vs dp p 3(bc - b^c^) 

which goes to infinity as t — > which is unsurprising given the divergence of the pressure 
at the sudden singularity. 

3.2 Puiseux and logarithmic series 

We shall now discuss some other series expansions around the singularity. Consider first 
a Puiseux series of the form 

oo 

a(t) = a J t3 + tn h + tl/qb i + * 2/ " & 2 + • • • ] 

3=0 

as an asymptotic form for the Friedmann-Lemaitre equations instead of eq. (jSJ). This 
suggests an ansatz for the quasi-isotropic form of the general inhomogeneous metric of 
the form 

lap = dap + K/3 t + C a/3 t n + ll) a p t m + L0 a/ 3 f H , 

10 



with 

\<n<m<r<2. 

If we perform the same analysis we will find again that a a p and b a p are arbitrary and 

c a /3 = -a a /3, and if> a p = —a a /3, and u a p = —a a p, 

with c, if), oj again arbitrary. The pressure and the density will be 

nbc „ i mlnb _ , rbu , , 
16tt P = 16vrp + — t"" 1 + _r t™" 1 + _ + , (11) 



16np = _ 2n (n-l)c tn _ 2 

and 



i m(m-l)^> ^ m _ n m(m-l)o; f _ n 
n(n — l)c n(n — l)c 



,(12) 



m — n > r — n > 0. 

Therefore we have the same phenomenon as with the series that was described in the 
previous section. If we keep adding terms of the form t s with 1 < s < 2 we will getting 
one more arbitrary function so if we add an infinite number of terms we will have an 
infinite number of functions. The reason that they are irrelevant is that they do not 
alter the nature of the singularity in its neighborhood as can be seen from the density 
and pressure eqs. (TTTT)-(fl2l). In other words, all three forms 

7«/3 — a a/3 + b a /3 t + C a/ 3 t U H , 

la/3 = Cta/3 + b a p t + C a /3 t n + lf) a/3 t™ + U a/3 f H , 

lap = a a(3 + b a pt + c a p t n + ip afS t m ln(t) • • ■ 



give 



2n(n — l)c 2 
limp = pO) and I6n limp = — lim t 

t^o t^o t->o 3 



and hence the same behaviour on approach to the sudden singularity as t — > 0. Thus 
the number of relevant independent functions is 9 = 6 + 6+1— 4 which is the number 
required for a specification of the general solution. 
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4 Discussion 



If we had sought a power-series expansion of the metric with time dependence dictated 
by the form of the radiation-dominated Friedmann-Lemaitre universe on approach to an 
initial big bang curvature singularity as t — > then it would have taken the form of eq. 
© with 

lap = ta a p + t 2 b af3 + ... (13) 

with a a p, and b a p spatial functions [301 ED 132] However, at most 3 arbitrary spatial 
functions in the series expansion are left free in this solution. A general solution of this 
form, with a radiation equation of state requires 8 independently arbitrary functions of 
space. The reason for the non-generic form of ( |T3|) for the radiation fluid is its quasi- 
isotropic character. On approach to a general relativistic cosmological singularity the 
anisotropic shear and 3-curvatures dominate the dynamics and destroy the quasi isotropic 
character of a radiation-dominated universe. 

By contrast, our results are quite different. We find that on approach to a late- 
time sudden singularity, where the density, expansion rate, and metric remain finite, 
it is possible for a quasi-isotropic solution to be part of the general solution of the 
Einstein equations. This solution has no equation of state and is characterised by nine 
independently arbitrary spatial functions. On approach to a sudden singularity the 
mean scale factor and the shear approach constant values, the chaotically anisotropic 
degrees of freedom which dominate at an initial curvature singularity are frozen out, and 
anisotropic velocities grind to a halt because of the divergent pressure and inertia. The 
situation is far simpler than for a initial vacuum or p < p dominated fluid singularity. 
Comparable simplicity is achieved in the non-singular late-time approach to a quasi- 
isotropic de Sitter universe in the presence of a cosmological constant or a p = —p 
fluid, shown by Starobinsky [35]. We have shown that there is an analogous 'no hair' 
behaviour on approach to a late-time sudden singularity and the general structure of 
such a finite-time singularity admits a simple quasi-isotropic form. 
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